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Abstract: A molecule’s absorption in the atmospheric infrared (IR) window (IRW) is an 
indicator of its efficiency as a greenhouse gas. A model for estimating the absorption of a 
fluorinated molecule within the IRW was developed to assess its radiative impact. This model 
will be useful in comparing different hydrofluorocarbons and hydrofluoroethers contribution to 
global warming. The absorption of radiation by greenhouse gases, in particular hydrofluoroethers 
and hydrofluorocarbons, was investigated using ab initio quantum mechanical methods. Least 
squares regression techniques were used to create a model based on this data. The placement and 
number of fluorines in the molecule were found to affect the absorption in the IR window and 
were incorporated into the model. Several group increment models are discussed. An additive 
model based on one-carbon groups is found to work satisfactorily in predicting the ab initio 
calculated vibrational intensities. 
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Introduction 


Greenhouse gases (GHGs), such as chlorofluorocarbons (CFCs), perfluorocarbons (PFCs), 
and hydro fluorocarbons (HFCs) can contribute significantly to climate change . 1 Since the 
Montreal Protocol and the initial IPCC assessment, the Environmental Protection Agency (EPA) 
has investigated the use of HFCs and hydrofluoroethers (HFEs) as possible replacements for 
CFCs . 2 Although the industrial use of CFCs have been phased out, PFCs and HFCs are currently 
used heavily in the air conditioning, carpeting, drug, and the electronics industries . 3 ’ 4 In addition, 
they are used as tracers in atmospheric studies and have some medical applications . 5 ' 6 These 
greenhouse gases (HFCs, PFCs), which possess carbon-fluorine bonds, are known to absorb 
heavily in the atmospheric infrared (IR) window. The relative climate impact of these gases, as 
expressed by the global warming potential (GWP), is sometimes thousands of times that of CO2 
due to their inherent ability to absorb radiation within the Earth’s atmospheric IR window 
(IRW). The GWP is defined as the ratio of the time-integrated radiative forcing from the 
instantaneous release of one kilogram of a trace substance relative to that of one kilogram of a 
reference gas, which is usually taken to be CO2 . 1 7 The radiative forcing of a molecule has been 
shown to depend on its radiative efficiency and atmospheric abundance . 7 A large molecular 
radiative efficiency has been shown by Bera et al. to arise when a molecule exhibits a large 
infrared absorption intensity within the IR window due to large dipole derivatives of highly polar 
bonds . 8 Moreover, HFCs and especially PFCs, tend to possess long atmospheric lifetimes since 
carbon-fluorine bonds are generally not very reactive with hydroxyl radicals in the atmosphere. 
All of these factors lead to a large global warming potential. Finally, based on the analysis of the 
molecular origin of the strong absorption of radiation within the atmospheric IR window , 8 design 
strategies were devised to minimize an industrial molecule’s radiative efficiency . 9 
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Due to the relatively shorter atmospheric lifetimes, HFEs have been recommended as 
alternatives for PFCs and HFCs. 10,11 HFEs generally have an increased reactivity with hydroxyl 
radicals, which decreases their atmospheric lifetimes. 12 However, they absorb IR radiation very 
efficiently in the Earth’s atmospheric window wavelength range. 913 The carbon-oxygen 
stretches fall within the atmospheric IRW and have a large IR intensity due to the polarity of the 
C-0 bond. This in combination with the carbon-fluorine vibrational stretch, which also occurs 
within the atmospheric IRW, makes perfluoroethers (PFEs) and HFEs very effective greenhouse 

8,9,12,14-18 

gases. 

For fluorinated hydrocarbons and ethers, molecular structure greatly impacts its ability to 
absorb IR radiation. 9 It has been shown that by strategic distribution of fluorine atoms along the 
carbon framework, the infrared radiation absorption capacity can be reduced by a factor of 
two. 8 ' 9 Based on this knowledge, a group increment scheme for absorption intensities is proposed 
and investigated in the present study. The idea of a group increment scheme is analogous to 
Benson’s work with molecular thermodynamic properties. Benson approached the idea of 
describing thermodynamic properties by distinguishing different groups within a molecule and 
applying additivity rules. 19 With the different groups and their corresponding thennodynamic 
quantities, one can simply add the values of the different groups within a given molecule and 
estimate one of several thermodynamic properties. In a similar manner, PFCs, HFCs, PFEs, and 
HFEs can be characterized by different groups present in them, such as CF 3 , CF 2 , CHF 2 , CH 2 F, 
CH 2 , CHF-0 etc, in order to estimate the integrated infrared absorption intensity for a given 
molecule. Each group has a ‘base value’, which is determined by a specific model. Each groups’ 
effective contribution, i.e. ‘base value’ towards the total infrared absorption, is then estimated 
using a regression analysis of infrared intensities on a test set of over 350 fluorocarbons that 
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included hydrofluorocarbons, perfluorocarbons, hydrofluoroethers and perfluoroethers in this 
work. Due to the nature of groups present in the partially or fully fluorinated hydrocarbons as 
opposed to partially and fully fluorinated ethers, two sets of base values are estimated for similar 
groups, or in other words, we develop a model for fluorocarbons and a separate model for 
fluoroethers. These ‘base values’ are expected to greatly simplify the process of estimating 
infrared absorption intensities and thus molecular radiative efficiencies for PFCs, HFCs, PFEs, 
and HFEs. Also, while the exact definition of the Earth’s IR window has varied from one work 
to another, for purposes of the present study we define the window as 775 cm" 1 to 1425 cm" 1 in 
order to be inclusive. We realize that this is a relatively large defined IRW, but we note that the 
vibrational frequencies we use are hannonic, and inclusion of anhannonic contributions would 
lower the stretching modes, in particular, by 15-30 cm' 1 . 

Before proceeding, we point out that there has been some discussion in the literature recently 
concerning a molecule’s radiative efficiency (RE) versus its integrated IR absorption intensity 
within the Earth’s IRW. 20,21 Firstly, we note that the term RE has been used inconsistently in the 
literature, and in some cases RE is used in a fashion in which it is essentially the integrated IR 
intensity. However, the most commonly used definition for RE is related to the molecule’s 
integrated IR intensity, but also is multiplied by the Planck function (usually using a narrow band 
model - see Refs. 22 and 23) 2223 , and can also take into account IR absorption by other 
molecules in the atmosphere. Essentially, the RE defined in this way is part of the radiative 
transfer model used as part of a climate model. One concern, however, that we have with this 
approach is that it does not take into account the number of photons that can be processed per 
second by a single molecule (i.e., the model lacks details from molecular physics). In any case, 
we plan to explore this further in the future, but this is beyond the scope of the present study. 
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Hence for the present investigation we limit the group increment scheme to the integrated IR 
intensity within the IRW, knowing full well that one could adapt this method for REs (computed 
using the narrow band model Planck function, etc.) and expect similar results. 

In the next two sections the methods for obtaining the ‘base values’ for the groups and their 
use in obtaining total absorption intensities are explained. 

Theoretical Methods 

The molecular geometries were optimized using the second-order Moller-Plesset perturbation 
theory (MP2) along with a double zeta plus polarization one-particle basis set, denoted DZP++, 
which includes diffuse functions on hydrogen atoms as well as the heavy atoms. 24 " 26 The MP2 
method was chosen rather than a density functional theory (DFT) approach, for example B3LYP, 
because while the B3LYP method may fortuitously yield harmonic frequencies that are closer to 
the experimental fundamental vibrational frequency compared to MP2, it is our experience that 
MP2 is superior in predicting the IR intensities. 27 ' 28 This is explained by the well known fact that 
DFT methods generally do not describe the diffuse part of wave functions reliably, and electrical 
properties such as the dipole moment, dipole derivatives, polarizabilities, etc. all depend 
significantly on the diffuse part of the wave function. All ab initio electronic structure 
calculations were perfonned using the Q-Chem3.1 quantum chemistry program. 29 Harmonic 
vibrational frequencies and IR intensities were computed at the same level of theory 
(MP2/DZP++). The IR intensities were computed under the double harmonic approximation, 
meaning the harmonic oscillator approximation was utilized in the mechanical and electrical 
framework. With this approximation, the IR vibrational intensities are proportional to the square 
of the dipole derivatives. 30 From comparison with experimental data, the IR intensities computed 
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at the MP2/DZP++ level of theory are expected to be between 5-20% larger than experimental 
values. 31 For detennining whether this type of model can work for IR intensities, the 
MP2/DZP++ level of theory is adequate. Linear least squares regression models were 
detennined using Minitab 15 Statistical Software (2007). 32 R 2 values, residuals, and fitted 
intensity values were computed using Minitab 15. 

Statistical Background 

The residual is defined as the difference between the intensity calculated by ab initio methods 
and the intensity predicted by the model, 

= y, - (!) 


The fitted intensities matrix that minimizes the sum of squares is given by the following formula 
where X represents the matrix that has n molecules as rows and m columns of x-variables. ft 
denotes the matrix consisting of the fitted parameters, and y represents the column matrix 
consisting of all ab initio calculated intensities. 

P-(X T X)~ l X T y (2) 

y = Xp = X(X T X)- 1 X T y (3) 

The coefficient of determination for a model with n molecules is used to assess the fit of the 
model. It is given in Eq. 4 below where y represents the mean of the ab initio calculated 
intensities. 


e = residual 





( 4 ) 
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y = ab initio calculated IR absorption intensity 

y = IR absorption intensity predicted by model 

X = matrix consisting of all input data for molecules 

= matrix consisting of least squares fitted parameters 
R 2 = coefficient of determination 

A summary of the linear regression analysis that we have used is contained in Ref. 33 . 33 

Results and Discussions 

B. Hydrofluorocarbon and Perfluorocarbon Models 

We have used ab initio quantum chemical methods to calculate the harmonic vibrational 
frequencies and intensities for IR absorption for partially and fully fluorinated two carbon, three 
carbon, four carbon molecules and partially fluorinated five-carbon molecules. For the two, 
three, and four carbon molecules, we started with singly fluorinated species and included all 
possible arrangements through molecules saturated with fluorine. From the ab initio computed 
intensities several model hydrocarbon and perfluorocarbons, shown in Table 1, were created. 

The hydrofluorocarbon and perfluorocarbon molecules are broken down into one-carbon 
groups. In the simplest models we assume that the base values for one-carbon groups are 
independent of their neighboring groups. The first two additive models for predicting the 
molecule’s absorption within the atmospheric IRW are based on this approximation and so do 
not include interaction terms between groups. Model 1 also includes an intercept term to account 
for any residual absorption within the atmospheric IRW, which is not accounted for by the group 
base values. When this regression was run, the CH 3 group was automatically removed from the 
regression due to high correlation. A group is considered highly correlated with other predictor 
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variables if the regression of that variable on the other predictor variables results in a regression 
equation where 1- R 2 is less than 10" 18 . In Model 2, no intercept term is allowed to enter and the 
integrated intensity is represented by the sum of base values of the constituent one-carbon 
groups. When Model 2 was run, the CH 3 group was no longer highly correlated with the other 
columns in the basis matrix and so was included in the final model. In a previous study by Bera 
et ah, it was shown that increasing the number of carbon-fluorine bonds leads to a nonlinear 
increase in IR intensity. 8 Since the carbon-hydrogen bond is less polar than the carbon-fluorine 
bond, the CH 3 group may have been highly correlated with other predictor variables because the 
contributions from the carbon-fluorine bonds overwhelm the contributions from a group that 
consists of only carbon-hydrogen bonds. However, the CH 2 group was not highly correlated with 
other predictor variables and was still included in the regression equation. Although it is 
included, the base values for the CH 2 group is much smaller than the base values for groups 
containing at least one carbon-fluorine bond. For both Model 1 and Model 2, each group is 
defined by the combination of atoms attached to a carbon. For example, the molecule 
CFH 2 CF 2 CFHCH 3 contains one CFH 2 group, one CF 2 group, one CFH group, and one CH 3 
group. The intensity in the IRW for this molecule was calculated by ab initio methods to be 
536.2 km mol". Models 1 and 2 are used in the following manner to predict this intensity: 

Model 1 

Intensity in IRW = 86.6 + 704 CF 3 + 303 CF 2 H + 77.2 CFH 2 + 289 CF 2 + 75.6 CFH + 14.8 CH 2 
= 86.6 + 704 (0) + 303 (0) + 77.2 (1)+ 289 (1) + 75.6 (1) + 14.8 (0) 

= 528.4 km mol" 1 

Model 2 
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Intensity in IRW = 748 CF 3 + 346 CF 2 H + 121 CFH 2 + 43.3 CH 3 + 289 CF 2 + 75.6 CFH + 14.8 
CH 2 

= 748 (0) + 346 (0) + 121 (1) + 43.3 (1) + 289 (1) + 75.6 (1) + 14.8 (0) 

= 528.9 km mol' 1 

Both Model 1 and 2 are in agreement with each other and predict the intensity in the 
atmospheric IRW well. In Model 1 and 2, the intensity in the IRW was under-predicted for the 
molecule CF 3 CH 2 CH 2 CF 3 . The histograms, Figure 1 and 2, of residuals for Models 1 and 2 are 
almost symmetric and bell-shaped, and show evidence of an outlier. For certain molecules that 
contain the same groups, Models 1 and 2 may predict the same intensities in IRW for different 
molecules. For example, with these models the molecules CF 3 CFHCH 2 CH 3 and 
CF 3 CH 2 CFHCH 3 would have the same predicted intensities. 

Models 3 and 4, included in the supporting information available through the website, are 
used in a similar way to Models 1 and 2. However, these two models contain terms that account 
for all possible interactions of groups within the molecule. In other words every one-carbon 
group feels the presence of all the other groups in the molecules. In Model 3, the CH 3 group and 
all CH 3 interaction terms were highly correlated with other predictor variables and were not 
included in the final model. In Model 4, without the constant tenn, all interaction terms involving 
CH 3 were highly correlated with other predictor variables and were automatically removed from 
the regression. Models 3 and 4 predict the same intensities in the IRW for a molecule. The 
interaction terms between the one-carbon groups are much smaller in magnitude than the one- 
carbon group values, and further, interaction tenns for groups separated by one carbon or more 
are even smaller. The histograms (included in the supporting material) of residuals for these 
models are slightly skewed left, but the highest abundance of residuals is near zero. In certain 
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cases, these models will predict the same intensities in IRW for different molecules. For 
instance, the molecules CF 3 CF 2 CH 2 CH 3 and CF 3 CH 2 CF 2 CH 3 would be predicted to have the 
same intensities in IRW when in reality their intensities differ. These models are used in the 
following way for the molecule CF 3 CH 2 CH 2 CF 2 H to predict the intensity of absorption within 
the IRW. This molecule has one CF 3 group, two CH 2 groups, one CF 2 H group, two CF 3 /CH 2 
interactions, two CF 2 H/CH 2 interactions, one CF 3 /CF 2 H interaction, and one CH 2 /CH 2 
interaction. 

Model 3 

Intensity in IRW = 99.5 + 680 CF 3 + 307 CF 2 H + 85.0 CFH 2 + 324 CF 2 + 63.2 CFH 

- 75.5 CH 2 - 12.9 CF 3 / CF 2 + 4.57 CF 3 / CFH + 60.0 CF 3 / CH 2 

+ 22.5 CF 3 / CF 3 + 3.9 CF 3 / CF 2 H - 3.4 CF 3 / CFH 2 - 20.7 CF 2 H / CF 2 

- 14.9 CF 2 H / CFH + 28.3 CF 2 H / CH 2 + 24. 1 CF 2 H / CF 2 H 
+ 2.3 CF 2 H / CFH 2 - 19.9 CFH 2 / CF 2 - 2.67 CFH 2 / CFH 

+ 22.8 CFH 2 / CH 2 - 6.9 CFH 2 / CFH 2 - 42.7 CF 2 / CF 2 - 9.5 CF 2 / CFH + 31.9 CF 2 / 
CH 2 + 22.2 CFH/CFH + 42.6 CFH/ CH 2 + 33.4 CH 2 / CH 2 

=99.5 + 680 (1) + 307 (1) + 85.0 (0) + 324 (0) + 63.2 (0) 

- 75.5 (2) - 12.9 (0) + 4.57 (0) + 60.0 (2) 

+ 22.5 (0) + 3.9 (1) - 3.4 (0) - 20.7 (0) 

- 14.9 (0) + 28.3 (2) + 24.1 (0) 

+ 2.3 (0)- 19.9(0) -2.67 (0) 

+ 22.8 (0) - 6.9 (0) - 42.7 (0) - 9.5 (0) + 31.9 (0) + 22.2 (0) + 42.6 (0) + 33.4 (1) 
= 1149.4 km mol " 1 
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Model 4 


Intensity in IRW = 730 CF 3 + 356 CF 2 H + 135 CFH 2 + 49.7 CH 3 + 324 CF 2 + 63.2 CFH 

- 75.5 CH 2 - 12.9 CF 3 / CF 2 + 4.57 CF 3 / CFH + 60.0 CF 3 / CH 2 

+ 22.5 CF 3 / CF 3 + 3.9 CF 3 / CF 2 H - 3.4 CF 3 / CFH 2 - 20.7 CF 2 H / CF 2 

- 14.9 CF 2 H / CFH + 28.3 CF 2 H / CH 2 + 24. 1 CF 2 H / CF 2 H 
+ 2.3 CF 2 H / CFH 2 - 19.9 CFH 2 / CF 2 - 2.67 CFH 2 / CFH 

+ 22.8 CFH 2 / CH 2 - 6.9 CFH 2 / CFH 2 - 42.7 CF 2 / CF 2 - 9.5 CF 2 / CFH 
+ 31.9 CF 2 / CH 2 + 22.2 CFH / CFH + 42.6 CFH / CH 2 + 33.4 CH 2 / CH 2 

= 730 (1) + 356 (1) + 135 (0) + 49.7 (0) + 324 (0) + 63.2 (0) 

- 75.5 (2) - 12.9 (0) + 4.57 (0) + 60.0 (2) 

+ 22.5 (0) + 3.9 (1) - 3.4 (0) - 20.7 (0) 

- 14.9 (0) + 28.3 (2) + 24.1 (0) 

+ 2.3 (0)- 19.9(0) -2.67 (0) 

+ 22.8 (0) - 6.9 (0) - 42.7 (0) - 9.5 (0) 

+ 31.9 (0) + 22.2 (0) + 42.6 (0) + 33.4 (1) 

= 1148.9 km moF 1 


Notice that in both of these models there are negative base values for some group interaction 
terms. This is due to the effects of multicollinearity, meaning that some of the predictor variables 
are highly correlated with the other predictor variables. For the CH 2 group, this does not mean 
that incorporating this group into the design of a molecule will decrease the absorption intensity 
in the IRW; it means some groups have a negative sign to balance the effects of high correlations 
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between predictor variables. Presence of the negative coefficients in these two models is 
undesirable. Inclusion of the interaction terms between the one-carbon groups has little effect on 
the usefulness of the models. 

Finally, since the interactions between groups that are separated by at least one carbon are 
weaker than the nearest neighbor interactions, Models 5 and 6 account for only the groups within 
the molecule and the nearest neighbor interactions. Model 5 has a constant tenn while Model 6 
does not. Both the models predict the same values for some molecules and are presented in SI. 
With these models, the molecule CF 2 HCFHCH 3 (which contains one CF 2 H group, one CFH 
group, one CH 3 group, one CF 2 H/CFH interaction, and one CH 3 /CFH interaction and by ab initio 
methods has an intensity in the IRW of 492.9 km mol" 1 ) would be entered in the following way: 
Model 5 

Intensity in IRW = - 14.9 + 760 CF 3 + 345 CF 2 H + 1 12 CFH 2 + 318 CF 2 + 1 12 CFH 
+ 24.0 CH 2 - 15.4 CF 3 / CF 2 - 42.9 CF 3 / CFH + 2.4 CF 2 H / CF 2 

- 25.6 CF 2 H / CFH + 1.6 CFH 2 / CF 2 - 0.8 CFH 2 / CFH + 58.1 CH 3 / CF 2 
+ 29.3 CH 3 /CFH - 61.5 CF 2 / CF 2 - 38.4 CF 2 / CFH - 14.9 CFH/ CFH 

= - 14.9 + 760 (0) + 345 (1)+ 112 (0) + 318 (0)+ 112 (1) 

+ 24.0 (0) - 15.4 (0) - 42.9 (0) + 2.4 (0) 

- 25.6 (1) + 1.6 (0) - 0.8 (0) + 58.1 (0) 

+ 29.3 (1) - 61.5 (0) - 38.4 (0) - 14.9 (0) 

= 445.8 km mol " 1 


Model 6 
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Intensity in IRW = 752 CF 3 + 338 CF 2 H + 104 CFH 2 - 7.5 CH 3 + 318 CF 2 + 1 12 CFH 
+ 24.0 CH 2 - 15.4 CF 3 / CF 2 - 42.9 CF 3 / CFH + 2.4 CF 2 H/ CF 2 

- 25.6 CF 2 H/ CFH + 1.6 CFH 2 / CF 2 - 0.8 CFH 2 /CFH + 58.1 CH 3 / CF 2 
+ 29.3 CH 3 / CFH - 61.5 CF 2 / CF 2 - 38.4 CF 2 /CFH - 14.9 CFH/CFH 

= 752 (0) + 338 (1)+ 104 (0) - 7.5 (1) + 318 (0) + 1 12 (1) 

+ 24.0 (0) - 15.4 (0) - 42.9 (0) + 2.4 (0) 

- 25.6 (1) + 1.6 (0) - 0.8 (0) + 58.1 (0) 

+ 29.3 (1) - 61.5 (0) - 38.4 (0) - 14.9 (0) 

= 446.2 km mol" 1 

Both models under predict the intensity in the IRW for the molecule CF 3 CH 2 CH 2 CF 3 . 
Disregarding outliers, the histograms (included in the supporting information) for these models 
are symmetric and bell-shaped curves centered on zero. These models also have some negative 
coefficients, which are once again due to the effects of multicollinearity. Despite being more 
sophisticated than the Models 1 and 2, these two models improve the prediction only marginally. 
Additionally the presence of the negative base values for certain groups makes them 
unappealing. 

In Models 1-6, the coefficient of determination is high, indicating that these models for 
hydrofluorocarbons and perfluorocarbons describe the data well. For terminal carbon groups in 
all of the models, the coefficient decreases as the number of carbon-fluorine bonds decreases. 
This is also true for the inner carbon groups. The base values trend indicates that placing more 
fluorine atoms on a carbon results in a nonlinear increase in infrared absorption intensity and 
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radiative efficiency. Of all the models, the Model 1 and 2, presented in Table 1, are the simplest 
and have the desirable quality of all positive base values, besides being accurate. Model 2, with 
no intercept, no interaction tenns and only base values for one-carbon units would be our method 
of choice for estimating the infrared intensities of hydro fluoro- and perfluorocarbons. 

C. Hydrofluoroether and Perfluoroether Models 

The second class of molecules we modeled are the hydrofluoroethers and perfluoroethers. 
We used ab initio quantum chemical methods to calculate the hannonic frequencies and 
vibrational intensities of IR absorption for all combinations of two carbon, three carbon, and four 
carbon ethers containing one fluorine atom through saturation with multiple fluorine atoms. 
From this data, the regressions in Table 2 were generated. 

Models 7 and 8 are used in a similar manner to models 1 and 2. However, these 
molecules include groups that contain oxygen atoms. In Models 7 and 8 , each one-carbon group 
is unaffected by the presence of neighboring groups. In the ether molecules, if a carbon group is 
next to an oxygen, that group is only counted once as a group with an oxygen. This means that a 
molecule such as CF 3 OCFHCFHCF 3 contains the four distinct groups CF 3 0, CFHO, CFH, and 
CF 3 . That is, the intensity in the atmospheric IRW would be predicted in the following way: 
Model 7 

Intensity in IRW = 132 + 838 CF 3 + 453 CF 2 H + 213 CFH 2 + 129 CH 3 + 301 CF 2 + 75.7 CFH + 
15.1 CH 2 + 1310 CF 3 0 + 752 CF 2 HO + 319 CFH 2 0 + 617 CF 2 0 + 251 CFHO 
= 132 + 838 (1) + 453 (0) + 213 (0) + 129 (0) + 301 (0) + 75.7 (1) + 15.1 (0) 

+ 1310 (1) + 752 (0) + 319 (0) + 617 (0) + 251 (1) 

=2606.7 km mol ' 1 
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Model 8 


Intensity in IRW = 904 CF 3 + 519 CF 2 H + 279 CFH 2 + 195 CH 3 + 301 CF 2 + 75.7 CFH +15.1 
CH 2 + 1376 CF 3 0 + 818 CF 2 HO + 385 CFH 2 0 + 65.9 CH 3 0 + 617 CF 2 0 + 251 
CFHO 

= 904 (1) + 519 (0) + 279 (0) + 195 (0) + 301 (0) + 75.7 (1) + 15.1 (0) + 1376 (1) 

+ 818 (0) + 385 (0) + 65.9 (0) + 617 (0) + 251 (1) 

=2606.7 km mol' 1 

In Model 7, the CH 2 0 group and the CH 3 0 were highly correlated with other predictor 
variables and were automatically removed from the regression equation. In Model 8, only the 
CH 2 0 group had high correlation with other predictor variables and was not included in the 
regression equation. In some extreme cases both Model 7 and Model 8 may predict the same 
intensities for two molecules, which are constitutional isomers with similar groups. For example, 
both Model 7 and Model 8 would predict the same intensities for CF 3 CF 2 OCH 2 CH 3 and 
CH 3 CF 2 OCH 2 CF 3 even though these molecules have slightly different absorption intensities and 
radiative efficiencies. 

Models 9 and 10, included in (supporting infonnation) are used for fluoroethers in a similar 
manner as Models 3 and 4 for fluorocarbons and are the most sophisticated models. They also 
contain oxygen groups, which are used in the same way as in Models 7 and 8, i.e. if a group is 
next to an oxygen atom then that group is only counted along with the oxygen. This model did 
not contain the possible interactions with terminal oxygenated groups since it was detennined 
from early models that interaction terms in general had fairly high correlations with other 
predictor variables and so not all interaction terms may be necessary. While the interaction terms 
do contribute to improving the additive model, it was decided that only some of the interaction 
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terms needed to be included. In Model 10 , when the regression was run the groups CH 2 0 , 
CH3/CF2, CH3/CFH, CH3/CH2, CF 2 0 /CH 2 , CFHO/CH2, CF3/CH2, CF2H/CH2, CFH2/CH2, 
CH2O/CH3, CH2O/CF2, CH2O/CFH, and CH2O/CH2 were highly correlated with other predictor 
variables and were removed from the regression equation. For Model 9, the CH3O group was 
highly correlated with other predictor variables in addition to the highly correlated terms from 
Model 10 and so all of these terms were removed from the regression. These two models also 
have the undesirable property of having negative interaction tenns between one-carbon groups. 

Model 1 1 is used in the same way as Model 7 except the constant term is multiplied by the 
number of carbons the molecule contains. In certain cases, Model 1 1 does predict the same 
intensities for two different molecules that contain the same groups. Models 7, 8, and 1 1 predict 
the same fitted intensities in the IRW for each molecule. Similarly, Model 12 is used like Model 
9 with the constant tenn multiplied by the number of carbons in the molecule. Models 9, 10, and 
12 give the same predicted intensities in the IRW for a molecule. All of the models have 
symmetric, bell shaped histograms that are centered on zero and have similar ranges for the 
residuals presented in the supporting information. The histograms in addition to the R 2 values 
indicate that Models 7-12 predict the intensities well for the hydro fluoro- and the 
perfluoroethers. The nearest-neighbor models do have a slightly lower average residual and do 
not predict the same intensities for different molecules, but are more complicated than the 
additive models. In addition, in Models 9-12 there are some negative coefficients. This is once 
again due to the effects of multicollinearity, since there are some predictor variables that are 
highly correlated with other predictor variables. The negative coefficients do not necessarily 
mean that the group or interaction does not positively contribute to the absorption intensity in the 
IRW, it means that in order to counteract the overestimation of effects of other dependent groups 
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some of the coefficients must be negative. In all of these models, the coefficients increase as 
more fluorines are added to a terminal or inner carbon. This is in agreement with previous 
studies, 8 ' 9 which show that an increase in carbon-fluorine bonds results in a nonlinear increase in 
intensity for absorption within the IRW, and also with the observations of Blowers et al. 15,16 
Also, the coefficients for the groups containing oxygen are higher than the coefficients of the 
equivalent non-oxygenated groups. This agrees with observations made previously that 
fluorinated carbons situated next to an oxygen result in a larger absorption in the IRW or a larger 
radiative efficiency. 9 ' 12 The additive models, especially Model 8, are most appealing for being 
simplest and having all positive base values and therefore will be favored by us for estimating 
the intensities of hydro fluoroethers and perfluoroethers. 

Conclusions 

Several models were created to assess the vibrational absorption intensity of a molecule’s 
absorption within the IRW. The additive models (models 1 and 2 for the hydrofluorocarbons and 
7 and 8 for the perfluorocarbons) gave good agreement with experiment while retaining positive 
coefficients. These are the models that we would recommend. The additive models (1, 2, 7 and 
8) account for the contributions, ‘base values’, from different one-carbon groups while the 
nearest- neighbor models (5, 6, 11, and 12) take into consideration the effects of interactions for 
groups that are next to each other. The interaction models (3, 4, 9, and 10) for 
hydrofluorocarbons and perfluorocarbons take into account all possible interactions between 
groups. The constant models and carbon constant models include a term that can account for 
residual intensities. The carbon constant and no-constant models include one term that was 
removed from the respective constant model due to a high correlation. The coefficients of the 
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models for both fluorinated ethers and fluorinated hydrocarbons indicate that increasing the 
number of fluorines on a carbon results in a nonlinear increase in the molecule’s absorption of 
radiation within the IRW. The ether models’ base values show that fluorinated carbons next to an 
oxygen atom contribute more to the IRW absorption intensity than their respective non- 
fluorinated counterparts do. The nearest neighbor models and interaction models include 
unphysical negative coefficients or base values; these negative coefficients are the result of fairly 
high correlations between predictor variables. All the models have high R 2 values indicating that 
they explain the data well. However, some of the models do predict the same intensities for 
different molecules. Overall, these models offer several simple ways to estimate the relative 
infrared absorption intensities and hence radiative efficiencies of fluorcarbons and fluoroethers. 

Finally, we note that use of these models together with climate models would necessitate that 
the climate model treat trace species that absorb radiation within the IRW separately and 
differently relative to more abundant atmospheric molecules (such as CO 2 , H 2 0, N 2 0, and CH 4 ). 
This would seem to be prudent since within the IRW there is, by definition, an over abundance 
of available photons to be absorbed, and that trace species will not be able to absorb and process 
all the photons of a given wavelength within the IRW. 
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Tables 


Table 1. Least Squares Regression Models for Hydrofluorocarbons and Perfluorocarbons. 
Averaged absolute residuals are expressed in km mol' 1 . 


Regression Model 

Constant 

terms? 

Quadratic 

Terms? 

Data 

Excluded 

R 2 

% 

Outliers 

Averaged 

Absolute 

Residuals 

Largest 
Residual 
km mol' 1 

Model 1 

Additive- With Intercept: 

Intensity in IRW = 86.6 + 
704 CF 3 + 303 CF 2 H + 
77.2 CFH 2 + 289 CF 2 + 
75.6 CFH + 14.8 CH 2 

Yes 

no 

All one- 

carbon 

molecules 

99.3 

cf 3 ch 2 ch 2 cf 3 

28.3 

138.4 for 

cf 3 ch 2 ch 2 cf 3 

Model 2 

Additive-No Intercept: 

Intensity in IRW = 748 
CF 3 + 346 CF 2 H+ 121 
CFH 2 + 43.3 CH 3 + 289 
CF 2 + 75.6 CFH + 14.8 

ch 2 

No 

no 

All one- 

carbon 

molecules 

99.3 

cf 3 ch 2 ch 2 cf 3 

28.3 

138.4 for 

cf 3 ch 2 ch 2 cf 3 


Figure 1. Histograms of Residuals for Models 1 and 2. (1) Histogram of residuals for Model 
1. (2) Histogram of residuals for Model 2. 
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Histogram of Residuals for Intensity in IRW With Intercept 



(l) 


( 2 ) 


Histogram of Residuals for Intensity in IRW Without Intercept 



Table 2. Least Squares Regession Models for Hydrofluoroethers and Perfluoroethers. 
Averaged absolute residuals are expressed in km mol" 1 . 


Regression Model 

Constant 

term? 

Quadratic 

Terms? 

R 2 

% 

Outliers 

Averaged 
Absolute 
Value of 
Residuals 

Largest 
Residual 
km mol' 1 

Model 7 

Additive W ith Constant: 

Intensity in IRW = 132 + 838 CF 3 + 453 
CF 2 H + 213 CFH 2 + 129 CH 3 + 301 CF 2 

+ 75.7 CFH+ 15.1 CH 2 + 1310 
CF 3 0 + 752 CF 2 HO + 319 CFH 2 0 

+ 617 CF 2 0 + 251 CFHO 

Yes 

no 

99.7 

Possible 

Outlier: 

CF 2 HOCF 2 CF 2 

cfh 2 

27.2 

113.3 for 
CF 2 HOCF 2 CF 2 

cfh 2 

Model 8 

Additive-No constant 

Intensity in IRW = 904 CF 3 + 5 19 CF 2 H + 

279 CFH 2 + 195 CH 3 + 301 CF 2 + 75.7 CFH 

+ 15.1 CH 2 + 1376 CF 3 0 + 818CF 2 H0 + 

385 CFH 2 0 + 65.9 CH 3 0 + 617 CF 2 0 + 251 

CFHO 

No 

no 

99.7 

Possible: 

CF 2 HOCF 2 CF 2 

cfh 2 

27.2 

113.3 for 
CF 2 HOCF 2 CF 2 

cfh 2 
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Figure 2. Histograms of Residuals for Hydrofluoroether and Perfluoroether Models. (7) 


Histogram of residuals for Model 7. (8) Histogram of residuals for Model 8. 
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